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Abstract 

The MSSM with large tan j3 and heavy squarks (Mg > 1 TeV) is a theoretically well mo- 
tivated and phenomenologically interesting extension of the SM. This scenario naturally 
satisfies all the electroweak precision constraints and, in the case of not too heavy slepton 
sector (Mg ^ 0.5 TeV), can also easily accommodate the (g — 2) M anomaly. Within this 
framework non-standard effects could possibly be detected in the near future in a few low- 
energy flavour- violating observables, such as B(B — > ru), B(B s d — > £ + £~), B(B — > X s j), 
and B(fi — > erf). Interpreting the (g — 2) M anomaly as the first hint of this scenario, 
we analyse the correlations of these low-energy observables under the additional assump- 
tion that the relic density of a Bino-like LSP accommodates the observed dark matter 
distribution. 



1 Introduction 

Within the Minimal Supersymmetric extension of the Standard Model (MSSM), the sce- 
nario with large tan/3 and heavy squarks is a particularly interesting subset of the pa- 
rameter space. On the one hand, values of tan /3 ~ 30-50 can allow the unification of top 
and bottom Yukawa couplings, as predicted in well- motivated grand-unified models [1]. 
On the other hand, heavy soft-breaking terms in the squark sector (both bilinear and tri- 
linear couplings) with large tan/3 and a Minimal Flavour Violating (MFV) structure [2,3] 
lead to interesting phenomenological virtues. On the one hand, this scenario can eas- 
ily accommodate all the existing constraints from electroweak precision tests and flavour 
physics. In particular, in a wide region of the parameter space, the lighetst Higgs boson 
mass is above the present exclusion bound. On the other hand, if the slepton sector is not 
too heavy, within this framework one can also find a natural description of the present 



(g — 2)^ anomaly. In the near future, additional low-energy signatures of this scenario 
could possibly show up in B[B — > ru), B(B sd — > and B(B — > X s ^) (see Ref. [4,5] 

for a recent phenomenological discussion). In the parameter region relevant to 5-physics 
and the (g — 2) M anomaly, also a few Lepton Flavor Violating (LFV) processes (especially 
fi — > ej) are generally predicted to be within the range of upcoming experiments. In 
this paper we analyse the correlations of the most interesting low-energy observables of 
this scenario under the additional assumption that the relic density of a Bino-like lightest 
supersymmetric particle (LSP) accommodates the observed dark matter distribution (the 
constraints and reference ranges for the low-energy observables considered in this work 
can be found in Sect. 13. 2p . 

Recent astrophysical observations consolidate the hypothesis that the universe is full 
of dark matter localized in large clusters [6]. The cosmological density of this type of 
matter is determined with good accuracy 

0.079 < tt CDM h 2 < 0.119 at2aC.L., (1) 

suggesting that it is composed by stable and weakly- interactive massive particles (WIMPs) . 
As widely discussed in the literature (see e.g. Ref. [7] for recent reviews), in the MSSM 
with .R-parity conservation a perfect candidate for such form of matter is the neutralino 
(when it turns out to be the LSP) [8]. In this scenario, due to the large amount of 
LSP produced in the early universe, the lightest neutralino must have a sufficiently large 
annihilation cross-section in order to satisfy the upper bound on the relic abundance. 

If the fi term is sufficiently large (i.e. in the regime where the interesting Higgs- 
mediated effects in flavour physics are not suppressed) and Mi is the lightest gaugino 
mass (as expected in a GUT framework), the lightest neutralino is mostly a Bino. Due 
to the smallness of its couplings, a Bino-like LSP tends to have a very low annihilation 
cross sectionj^ However, as we will discuss in Section |2j in the regime with large tan (3 
and heavy squarks the relic-density constraints can easily be satisfied. In particular, the 
largest region of the parameter space yielding the correct LSP abundance is the so-called 
A-funnel region [9]. Here the dominant neutralino annihilation amplitude is the Higgs- 
mediated diagram in Fig. [TJ Interestingly enough, in this case several of the parameters 
which control the amount of relic abundance, such as tan (3 and the heavy Higgs masses, 
also play a key role in flavour observables. As a result, in this scenario imposing the 
dark-matter constraints leads to a well-defined pattern of constraints and correlations on 
the low-energy observables which could possibly be tested in the near future. The main 
purpose of this article is the investigation of this scenario. 

The interplay of (g — 2) M , B(B S ^ — > B[B — > X s j), and dark-matter constraints 

in the MSSM have been addressed in a series of recent works, focusing both on relic 
abundance [10] and on direct WIMPs searches [11]. Our analysis is complementary to 
those studies for two main reasons: i) the inclusion of BiB — > tu), which starts to play a 
significant role in the large tan f3 regime, and will become even more significant in the near 

1 If the conditions on \i and M\ are relaxed, the LSP can have a dominant Wino or Higgsino component 
and a naturally larger annihilation cross-section. This scenario, which is less interesting for flavour physics, 
will not be analysed in this work. 
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Figure 1: Higgs- mediated neutralino annihilation amplitude. 



future; ii) the study of a phenomenologically interesting region of the MSSM parameter 
space which goes beyond the scenarios analysed in most previous studies (see Section [2]). 

The plan of the paper is the following: in Section [2] we recall the ingredients to evaluate 
the relic density in the MSSM, and determine the key parameters of the interesting A 
funnel region. In Section[3]we present a brief updated on the low-energy constraints on this 
scenario; we analyse constraints and correlations on the various low-energy observables 
after imposing the dark-matter constraints; we finally study the possible correlations 
between (g — 2) M and the lepton-flavour violating decays B{ji — > ej) and B(t — > fij). The 
results are summarized in the Conclusions. 



2 Relic Density 

In the following we assume that relic neutralinos represent a sizable fraction of the ob- 
served dark matter. In order to check if a specific choice of the MSSM parameters is 
consistent with this assumption, we need to ensure two main conditions: i) the LSP is 
a thermally produced neutralino; ii) its relic density is consistent with the astrophysical 
observation reported in Eq. ([I]). 

In the MSSM there are four neutralino mass eigenstates, resulting from the admixture 
of the two neutral gauginos (W°, B) and the two neutral higgsinos (H®, H®)- The lightest 
neutralino can be defined by its composition, 

Xi = Z U B + Z 12 W° + Z XZ H\ + Z U H° 2 (2) 

where the coefficients Zu and the mass eigenvalue (M^) are determined by the diagonal- 
ization of the mass matrix 

— rriz cos 3 Sw Tfiz sin (3 Sw \ 

M 2 m z cos (3 cw —m z sin (3 cw 

niz cos (3 cw — // 

-m z sin (3 cw —fJ> / 

As usual, $w denotes the weak mixing angle (cw = cos 6^, sw = sm 9w) and (3 is defined 
by the relation tan/? = t^/vi, where t>2(i) is the vacuum expectation value of the Higgs 
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— mz cos (3 sw 
V m z sin 3 s w 
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coupled to up(down)-type quarks; Mi and M 2 are the soft-breaking gaugino masses and 
\i is the supersymmetric-invariant mass term of the Higgs potential. 

In order to compute the present amount of neutralinos we assume a standard thermal 
history of the universe [12] and evaluate the annihilation and coannihilation cross-sections 
using the micrOMEGAs [13] code. Since we cannot exclude other relic contributions in 
addition to the neutralinos, we have analysed only the consistency with the upper limit in 
Eq. ([1]). This can be translated into a lower bound on the neutralino cross sections: the 
annihilation and coannihilation processes have to be effective enough to yield a sufficiently 
low neutralino density at present time. 

With respect to most of the existing analysis of dark-matter constraints in the MSSM, 
in this work we do not impose relations among the MSSM free parameters dictated by 
specific supersymmetry-breaking mechanisms. Consistently with the analysis of Ref. [4], 
we follow a bottom-up approach supplemented by few underlying hypothesis, such as 
the large value of tan/3 and the heavy soft-breaking terms in the squark sector. As 
far as the neutralino mass terms are concerned, we employ the following two additional 
hypotheses: the GUT relation M\ « M 2 /2 « M 3 /6, and the relation fi > M 1; which 
selects the parameter region with the most interesting Higgs-mediated effects in flavour 
physics (see Section [3]) H These two hypotheses imply that the lightest neutralino is Bino- 
like (i.e. Zyy Zi^i) with a possible large Higgsino fraction when \i = O(Mi). Due to the 
smallness of the B couplings, some enhancements of the annihilation and coannihilation 
processes are necessary in order to fulfill the relic density constraint. In general, these 
enhancements can be produced by the following three mechanisms [7, 14] : 

• Light sfermions. For light sfermions, the t-channel sfermion exchange leads to a 
sufficiently large annihilation amplitude into fermions with large hypercharge. 

• Coannihilation with other SUSY particles. If the next-to-lightest supersymmetric 
particle (NLSP) mass is closed to M^, the coannihilation process NLSP+LSP — > SM 
can be efficient enough to reduce the amount of neutralinos down to the allowed 
range. A relevant coannihilation process in our scenario occurs when the NLSP 
is the lightest stau lepton (stau annihilation region). This mechanism becomes 
relevant when the lightest stau mass, M~ R w M~ — m T fi tan /3, satisfies the following 
condition 

Ma < M fR < 1.1 x Mfc . (4) 

Other relevent coannihilation processes take place when \x is sufficiently close to 
Mi. In this case the LSP coannihilation with a light neutralino or chargino (mostly 
higgsino-like and thus with mass M-o ^+ ~ (/,), can become efficient. 

2 These two assumptions are not strictly necessary. From this point of view, our analysis should not be 
regarded as the most general analysis of dark-matter constraints in the MSSM at large tan j3. We employ 
these assumptions both to reduce the number of free parameters and to maximize the potentially visible 
non-standard effects in the flavour sector. In particular, the condition fi > Mi does not follow from 
model-building considerations (although well-motivated scenarios, such as mSUGRA, naturally predict 
/i > Mi in large portions of the parameter space), rather from the requirement of non- vanishing large- 
tan (3 effects in B — > and other low-energy observables [15-17] (which provide a distinctive signature 
of this scenario) . 
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Figure 2: Allowed regions in the M\-Mh plane satisfying the relic density constraint 
Vth 2 < 0.119 for Mq = 2M i = \A V \ = 1 TeV. Left panel: /i = 0.5 TeV with tan/5 = 20 
(green), 30 (green+red) and 50 (all points up to Mi « 480 GeV, see right panel). Right 
panel: /i = 0.5 TeV (black) and \i = 1 TeV (blue) for tan/? = 50 (a color version of this 
and the following figures can be found in the on-line version of this article). 



• Resonant processes. Neutralinos can efficiently annihilate into down-type fermion 
pairs through s-channel exchange close to resonance (see Fig.[T]). At large tan/?, the 
potentially dominant effect is through the heavy-Higgs exchange (A and H°) and 
in this case the resonant condition implies 

4% « M A /2 . (5) 

At resonance the amplitude is proportional to (M^/M^) (^11^13,14) (ra& (T /rajy) tan/? 
which shows that the lightest neutralino must have a non-negligible higgsino com- 
ponent (^13,14 7^ 0), and that the annihilation into b and r fermions grows at large 
tan/3 (relaxing the resonance condition). 

Because of the heavy squark masses, the first of these mechanisms is essentially ex- 
cluded in the scenario we are considering: we assume squark masses in the 1 — 2 TeV range 
and, in order to maintain a natural ratio between squark and slepton masses, this implies 
sleptons masses in the 0.3-1 TeV range. The second mechanism can occur, but only in 
specific regions. On the other hand, the s-channel annihilation \X ~~ * H, A — > bb{r + T~) 
can be very efficient in a wide region of the parameter space of our scenario. 

In Fig. [2]we explore the dark-matter constraints in the M\-Mn plane, assuming heavy 
squarks and sleptons {Mq = 1 TeV, = 0.5 TeV) and large trilinear couplings (\Ajj\ = 1 
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TeV). The allowed points have been obtained for different values of /i and tan/3. The 
dependence on tan/3 at fixed /i (/i = 0.5 TeV) is illustrated by the left panel, while 
the /i dependence at fixed tan/? (tan/3 = 50) is illustrated by the right panel. In all 
cases the heavy-Higgs resonant region, Mi « Mh,a/2, is the most important one @ The 
M#-independent regions for Mi > 450 GeV and Mi ~ 60 GeV are generated by the x 
coannihilation mechanisms and the h resonance amplitude, respectively. As can be seen, 
in the heavy-Higgs resonant case the allowed region becomes larger for larger M# values: 
this is because the Higgs width grows with M# and therefore the resonance region becomes 
larger. For a similar reason, and also because the annihilation cross sections grow with 
tan /3, the allowed region becomes larger for larger tan (5 values. As far as the \x dependence 
is concerned, the heavy-Higgs resonance region is larger for small /i values. This is because 
the xxA coupling, relevant in the resonant process, depends on the Higgsino component 
of x' f° r large /i, x is almost a pure B and the xxA coupling is suppressed. This fact can 
be used to set a theoretical upper limit on the \i parameter in this specific framework: jj, 
must be larger than Mi in order to reproduce a Bino LSP, but it should not be too heavy 
not to suppress too much the Bino annihilation amplitude. 

Notice that in the right panel of Fig. [2] only the x~ f coannihilation process is active 
when fj, = 1 TeV. On general grounds, given a left-handed slepton mass M~ tl the stau 
coannihilation region appears for lower Mi if /i increases, since Mf decreases with increas- 
ing /J,. Notice also that the h resonance region disappears for large /i, due to the smallness 
of the xxh coupling. In both figures points with Mf < M^ have not been plotted since 
they are ruled out. 

In summary, the MSSM scenario we are considering is mainly motivated by flavour- 
physics and electroweak precision observables. As we have shown in this Section, in this 
framework the dark matter constraints can be easily fulfilled with a Bino-like LSP and 
an efficient Higgs-mediated Bino annihilation amplitude. The latter condition implies a 
strong link between the gaugino and the Higgs sectors (most notably via the relation 
Mi w Mff/2). This link reduces the number of free parameters, enhancing the possible 
correlations among low-energy observables. 

3 Low-energy observables 

In this Section we analyse the correlations of new-physics effects in = (g — 2) M /2, 
B(B — > ru), B(B s ^ d — > £ + £~), B(B — > X/y), B(fM — > &y), and £>(r — > /ry), after imposing 
the dark matter constraints. As far as the I?-physics observables are concerned, we use 
the existing calculations of supersymmetric effects in the large tan/? regime which have 

3 We recall that for sufficiently heavy Mjj > 300GeV, the heavy Higgses are almost degenerate: 
Mh ~ Ma- We also recall that, within mSUGRA models, Mh, Mi and tan/3 are not independent 
parameters. In this case, the A-funnel condition Mh ~ 2 Mi is achieved only in the very large tan/? 
regime 45 < tan/3 < 60. In our scenario, where Mh and Mi are assumed to be free parameters, this 
constraints is relaxed and smaller values of tan/3 are also allowed. 
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been been recently reviewed in Ref. [4, 5]o However, since a few inputs have changed 
since then, most notably the BiB — > tv) measurements [22,23] and the SM calculation of 
BiB — > X s -f) [25], in the following we first present a brief updated on these two inputs. 
We then proceed analysing the implications on the MSSM parameter space of a M and 
.B-physics observables after imposing the dark matter constraints. Finally, the possible 
correlations between a M and the lepton-flavour violating decays £>(/! — > ej) and Bij — ► /ry) 
in this framework are discussed. 



3.1 Updated constraints from B — > rv and B —> X s j 

Due to its enhanced sensitivity to tree-level charged-Higgs exchange, [19] B — > rv is one 
of the most clean probes of the large tan f3 scenario. The recent 5-factory results [22,23] , 

B[B -> rz/) Babar = (0.88lg;gf(stat) ± O.ll(syst)) x 10~ 4 , 
B{B - ™) Belle = (1.79l^(stat)l^(syst)) x 10- 4 , (6) 

leads to the average BiB — ► rz/) cxp = (1.31 ± 0.49) x 10~ 4 . This should be compared 
with the SM expectation B(B -> rz/) SM = G|m B m^/i|K 6 | 2 (l-m^/m|) 2 /(87rr B ), whose 
numerical value suffers from sizable parametrical uncertainties induced by /# and V^,. 
According to the global fi|] of Ref. [24], the best estimate is B(B — > rz/) SM = (1.41 ± 
0.33) x 10 -4 , which implies 

B ex P(B — > rz/) 

^ = gSM( g ^ ry ) = 0-93 ±0-41. (7) 

A similar (more transparent) strategy to minimize the error on BiB — > ru) SM is the direct 
normalization of BiB — > rv) to AM# d , given that Bd-Bd is not affected by new physics in 
our scenario [4]. In this case, using BB d (rrib) = 0.836±0.068 and |V^/Vtd| = 0.473±0.024 
[24], we get 

,ex P _ B^jB -> TV)/AM% 

= 1.27 ±0.50 = 1.27±0.48 cxp ±0.10| BSd |±0.13|y WVtd | , (9) 

in reasonable agreement with Eq. (J7J). Although perfectly compatible with 1 (or with 
no new physics contributions), these results leave open the possibility of O(10% — 30%) 



4 See in particular Ref. [18] for B -> X sll Ref. [15-17] for B(B s>d -> Ref. [4,19] for B(B -» tu), 

and Ref. [20] for (g — 2)^/2. After this work was completed, a new theoretical analysis of large tan /3 effects 
in B physics, within the MFV-MSSM, has appeared [21]. As shown in Ref. [21], the renormalization of 
both tan (3 and the Higgs masses may lead to sizable modifications of the commonly adopted formulae 
for AMb s d (see Ref. [17]), which are valid only in the Mh ^> raw limit [3]. On the numerical side, 
these new effects turn out to be non-negligible only in a narrow region of light Mh {Ma ~ 160 GeV or 
Mh ~ 180 GeV) which is not allowed within our analysis. These new effects are therefore safely negligible 
for our purposes. 

In Ref. [24] the value of Jb is indirectly determined taking into account the information from both 
Bd~Bd and B s -B s mixing. 



7 



negative corrections induced by the charged-Higgs exchange. The present error on Rb tv 
is too large to provide a significant constraint in the MSSM parameter space. In order 
to illustrate the possible role of a more precise determination of B cxp (B — ► tu), in the 
following we will consider the impact of the reference range 0.8 < Rbtv < 0.9. In the 
next 2-3 years, at the end of the I?-factory programs, we can expect a reduction of the 
experimental error on B(B — ► tv) of a factor of 2-3. Depending on the possible shift of 
the central value of the measurement [note the large spread among the two central values 
in Eq. (JSJ)] the upper bound Rbtv < 0.9 could become the true 68% or 90% CL limit. 

The B — > Xg'-f transition is particularly sensitive to new physics. However, contrary to 
B — > tv, it does not receive tree-level contributions from the Higgs sector. The one-loop 
charged-Higgs amplitude, which increases the rate compared to the SM expectation, can 
be partially compensated by the chargino-squark amplitude even for squark masses of 
0(1 TeV). According to the recent NNLO analysis of Ref. [25], the SM prediction is 

B(B -> X sT , £ 7 > 1.6 GeV) SM = (3.15 ± 0.23) x 10~ 4 , (10) 
to be compared with the experimental average [26-28] 

B[B -> X sT , £ 7 > 1.6 GeV)) cxp = (3.55 ± 0.24) x 10" 4 . (11) 
Combining these results, we obtain the following la CL interval 

< R »*> ~ B^(B^X sl ) < L24 (12) 

which will be used to constrain the MSSM parameter space in the following numerical 
analysis^ 



3.2 Combined constraints in the MSSM parameter space 

The combined constraints from low-energy observables and dark matter in the tan /3—Mh 
plane are illustrated in Figures [3] and HI The plots shown in these figures have been 
obtained setting M- q = 1.5 TeV, \A V \ = 1 TeV, \i = 0.5 or 1 TeV, and M z = 0.4 or 
0.3 TeV. The two sets of figures differ because of the sign of A v . The gaugino masses, 
satisfying the GUT condition M2 ~ 1M\ « Ms/3, have been varied in each plot in order to 
fulfill the dark- matter conditions discussed in the previous Section (see Figure [2]). These 
conditions cannot be fulfilled in the gray (light-blue) areas with heavy Mjj, while the 
yellow band denotes the region where the stau coannihilation mechanism is active. The 
remaining bands correspond to the following constraints/reference-ranges from low-energy 
observables^ 

6 A slightly larger (and less standard) range is obtained taking into account the corrections associated 
to the E 1 cut in Ref. [29]. For simplicity, in our numerical analysis we have used Eq. (Tl2|) as reference 
range. The B —> X s ^ rate in the MSSM has been evaluated using the approximate numerical formula of 
Ref. [30], which partially takes into account NNLO effects. 

7 For the sake of clarity, the resonance condition Mjj = 2M\ has been strictly enforced in the bands 
corresponding to the low-energy observables. Similarly, the stau coannihilation region has been deter- 
mined imposing the relation 1 < Mf R /M^ < 1.1. 
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• B — > X s 7 [1.01 < Rss-y < 1-24]: allowed region between the two blue lines. 

• a M [2 < 10 9 (a^ xp — ) < 4 [31]]: allowed region between the two purple lines. 

• B — > /i 4 ' a~ [B exp < 8.0 x 10~ 8 [32]]: allowed region below the dark-green line. 

• AMb s [AMb 3 = 17.35 ± 0.25 ps" 1 [33]]: allowed region below the gray line. 

• B — > tv [0.8 < Rbtu < 0.9]: allowed region between the two black lines [ red (green) 
area if all the other conditions (but for a M ) are satisfied]. 

In the excluded regions at large Mh (light-blue areas) the neutralino cannot satisfy 
the resonance condition ~ Mff/2 and, at the same time, be lighter than the sleptons. 
This is why the excluded regions become larger for lighter M~ t For the same reason, 
the excluded regions become larger for larger values of fi (we recall that M~ R « M~ — 
m T fi tan (3). We stress that in all cases we have explicitly checked the consistency with 
electroweak precision tests and the compatibility with exclusion bounds on direct SUSY 
searches. By construction, these conditions turn out to be naturally satisfied in the 
scenarios we have considered. The mot delicate constraint is the value of the lightest 
Higgs boson mass (mh), which lies few GeV above its exlusion bound. In particular, we 
find 118 GeV < m h < 120 GeV in the plots of Figure El and 117 GeV < m h < 119 GeV 
in Figure HI 

As can be seen, in Figure [3] the B — > constraint is always easily satisfied for 
Mb > 300 GeV, or even lighter Mb values for large tan/? values. This is because the 
new range in Eq. (TT2"|) allows a significant (positive) non-standard contribution to the 
B — > Xg'j rate. Moreover, having chosen Ay < 0, the positive charged- Higgs contribution 
is partially compensated by the negative chargino-squarks amplitude. In Figure HI where 
Ajj > 0, the B — > Xg'-f constraints are much more stringent and almost tan /3-independent. 
It is worth noting that in Figure [3] the B — > Xgj information also exclude a region at 
large Mh'- this is where the chargino-squarks amplitude dominates over the charged-Higgs 
one, yielding a total negative corrections which is not favored by data. As already noted 
in [4], the precise AMb b measurement and the present limit on B — > do not pose 

any significant constraint. 

A part from the excluded region at large Mb, the most significant difference with 
respect to the analysis of Ref. [4] (where dark-matter constraints have been ignored) is 
the interplay between and 5-physics observables. The correlation between Mi and 
Mb imposed by the dark matter constraint is responsible for the rise with Mb of the 
a M bands in Figures [3] and HI This makes more difficult to intercept the B Xgj and 
B — > tv bands and, as a result, only a narrow area of the parameter space can fulfill all 
constraints. In particular, with the reference ranges we have chosen, the best overlap 
occurs for moderate/large values of tan/3 and low values of /x and M%. 

On the other hand, we recall that the B — ► tv band in Figure [3] does not correspond 
to the present experimental determination of this observable, but only to an exemplifying 
range. Assuming a stronger suppression of B(B — > tv) with respect to its SM value would 
allow a larger overlap between the B X s 7 and B — ► tv bands in the regions with higher 



9 



Mm 




200 400 600 ., 800 1000 1200 200 400 600 ., 800 1000 1200 

M H M H 



Figure 3: Combined constraints from low-energy observables and dark matter in the 
tan (3-Mh plane. The plots have been obtained for Mg = 1.5 TeV Ajj = — 1 TeV, and 
\fi,M$ = [1.0,0.4] TeV (upper left); \ji,Mg] = [0.5,0.4] TeV (upper right); [/i, M £ ~] = 
[1.0,0.3] TeV (lower left); [/i, M E ] = [0.5,0.3] TeV (lower right). The light-blue area is 
excluded by the dark-matter conditions. Within the red (green) area all the reference 
values of the low-energy observables (but for a M ) are satisfied. See main text for more 
details. The yellow band denote the area where the stau coannihilation mechanism is 
active (1 < M fR /M^ < 1.1); in this area the A-funnel region and the stau coannihilation 
region overlap. 
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Figure 5: Aa^ = (g^ — g^ M )/2 vs. the slepton mass within the funnel region taking 
into account the B — > X s j constraint and setting Rb T u > 0.7 (blue), Rbtu > 0.8 (red), 
Rbtv > 0.9 (green). The supersymmetric parameters have been varied in the following 
ranges: 200 GeV < M 2 < 1000 GeV, 500 GeV < fj, < 1000 GeV, 10 < tan/3 < 50. 
Moreover, we have set Ajj = — 1 TeV, Mq = 1.5 TeV, and imposed the GUT relation 
Mi « M 2 /2 « M3/6. 

values of tan /3, /i and M|. While if the B(B — ► rz/) measurement will converge toward 
the SM value, for the reference values of \i and M~ t chosen in the figures (fi > 0.5 GeV, 
Mg > 0.3 GeV) we deduce that: i) for Rb T u > 0.8 the non-standard contribution to a M 
cannot not exceed 3 x 10~ 9 ; ii) for Rb T v > 0.9 the non-standard contribution to a M cannot 
not exceed 2 x 10~ 9 . An illustration of how the non-standard contribution to a M varies as 
a funtion of M^, imposing different bounds on Rb T u, is shown in Figure [5j Moreover, if 
the B(B — > tv) measurement will converge toward the SM value and the a M constraint is 
not considered, the green areas in Figures [3] and @] are enlarged, allowing also lower tan (5 
values. 

In short, the main result of this analysis is that in a scenario with heavy squarks 
and large trilinear couplings, the constraints and reference ranges for the low-energy 
observables described above favor a charged Higgs mass in the 400 — 600 GeV range and 
tan f3 values in the 20-40 range. The structure of the favored tan f3 — M H region depends 
on other SUSY parameters, mainly fi and M|. Lower slepton masses shift the region 
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toward lower Mh and lower tan/? values (in order to reproduce the (g — 2) M anomaly and 
a neutralino LSP), while large \l values reduce the favored region selecting larger Mu and 
tan/5 values. 

The analysis of future phenomenological signals of this scenario at LHC and other 
experiments is beyond the scope of this work. However it should be noticed that both 
squarks and gluinos are rather heavy (around or above 1 TeV) and therefore not easily 
detectable. On the other hand, a direct detection of the charged Higgs and/or of the 
sleptons should be possible. In this case, the combination of high-energy and low-energy 
observables would allow to determine the tan j3 parameter very precisely. 

3.3 Correlation between LFV decays and [g — 2)^ 

As we have seen from the analysis of Figures [3] and HJ a key element which characterizes the 
scenario we are considering is the interplay between (g — 2) M and 5-physics observables. 
Since (g — 2) M is affected by irreducible theoretical uncertainties [31], it is desirable to 
identify additional observables sensitive to the same (or a very similar) combination of 
supersymmetric parameters. An interesting possibility is provided by the LFV transitions 
£i — > ij'-f and, in particular, by the /i — ► ey decay. Apart from the unknown overall 
normalization associated to the LFV couplings, the amplitude of these transitions are 
closely connected to those generating the non-standard contribution to a M [34]. 

LFV couplings naturally appear in the MSSM once we extend it to accommodate the 
non-vanishing neutrino masses and mixing angles by means of a supersymmetric seesaw 
mechanism [35]. In particular, the renormalization-group-induced LFV entries appearing 
in the left-handed slepton mass matrices have the following form [35]: 



where Y v are the neutrino Yukawa couplings and c v is a numerical coefficient, depending 
on the SUSY spectrum, typically of (9(0.1-1). As is well known, the information from 
neutrino masses is not sufficient to determine in a model-independent way all the seesaw 
parameters relevant to LFV rates and, in particular, the neutrino Yukawa couplings. To 
reduce the number of free parameters specific SUSY-GUT models and/or flavour sym- 
metries need to be employed. Two main roads are often considered in the literature (see 
e.g. Ref. [36] and references there in): the case where the charged- lepton LFV couplings 
are linked to the CKM matrix (the quark mixing matrix) and the case where they are 
connected to the PMNS matrix (the neutrino mixing matrix). These two possibilities 
can be formulated in terms of well-defined flavour-symmetry structures starting from the 
MFV hypothesis [37,38]. A useful reference scenario is provided by the so-called MLFV 
hypothesis [37], namely by the assumption that the flavour degeneracy in the lepton sector 
is broken only by the neutrino Yukawa couplings, in close analogy to the quark sector. 
According to this hypothesis, the LFV entries introduced in Eq. ([TBI assume the following 




(13) 
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form 



LL 



v)%3 



atm yr 
v 2 



(14) 



where M VR is the average right-handed neutrino mass and U denote the PMNS matrix. 

Once non-vanishing LFV entries in the slepton mass matrices are generated, LFV rare 
decays are naturally induced by one-loop diagrams with the exchange of gauginos and 
sleptons (gauge-mediated LFV amplitudes). H In particular, the leading contribution due 
to the exchange of charginos, leads to 



48vr 3 a 



a 2 UM % \ ho (Ml/Mf, n 2 /Mj 
(M|-^) 



47T 



ml 



5^tan/3, 



(15) 



where the loop function f2 C (x,y) is defined as f2c(x,y) = f2c(x) — f2c(y) in terms of 

-a 2 - 4a + 5 + 2(2a + 1) In a 



/2 C (a) 



2(1 



(16) 



Given that both ^ — > 1^ and Aa M = (g M — g^ M )/2 are generated by dipole operators, 
it is natural to establish a link between them. To this purpose, we recall the dominant 
contribution to Aa M is also provided by the chargino exchange and can be written as 



a a2 2 



^M 2 \92c [MllMl^lM'i 



VTLi 



tan/? . 



(17) 



with g C 2(x,y) defined as f C 2(x,y) in terms of 

(3 - 4a + a 2 + 2 log a) 



(18) 



It is then straightforward to deduce the relation 



B(£i - la) 



48tt 3 « 

Or p 



Aa, 



f 2c (Mj/Ml^/M? 
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°LL 



(19) 



To understand the relative size of the correlation, in the limit of degenerate SUSY spec- 
trum we get 



Aa, 



20 x 10 



-10 



x 



1 x 10- 4 |oii| 2 

23 |2 
LL\ 



2 x 10~ 5 \5 



T 



tA ■ 



(20) 



8 An additional and potentially large class of LFV contributions to rare decays comes from the Higgs 
sector through the effective LFV Yukawa interactions induced by non-holomorphic terms [43]. However, 
these effects become competitive with the gauge-mediated ones only if tan/3 ~ O(40 — 50) and if the 
Higgs masses are roughly one order of magnitude lighter then the slepton masses [44] . Since we consider 
a slepton mass spectrum well below the TeV scale, Higgs mediated LFV effects do not play a relevant 
role in our analysis. 
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Figure 6: Expectations for B{p, — > ej) and B(r — > fij) vs. Aa^ = (g^ — g^ M )/2, 
assuming \5 L 2 L \ = 10~ 4 and \S^ L \ = 10~ 2 . The plots have been obtained employ- 
ing the following ranges: 300 GeV < M~ e < 600 GeV, 200 GeV < M 2 < 1000 GeV, 
500 GeV < ii < 1000 GeV, 10 < tan/3 < 50, and setting A v = -1 TeV, M~ q = 1.5 TeV. 
Moreover, the GUT relations M 2 ~ 1M\ and M3 « 6M1 are assumed. The red areas cor- 
respond to points within the funnel region which satisfy the I?-physics constraints listed 
in Section [3J2] [B(B 8 //+//") < 8 x 10" s , 1.01 < R Bsy < 1.24, 0.8 < R Brv < 0.9, 
AM Bs = 17.35 ±0.25 ps -1 ]. 

A more detailed analysis of the stringent correlation between the — > ij^y transitions and 
Acin in our scenario is illustrated in Figj6j Since the loop functions for the two processes 
are not identical, the correlation is not exactly a line; however, it is clear that the two 
observables are closely connected. We stress that the numerical results shown in Figj6j 
have been obtained using the exact formulae reported in Ref. [39] for the supersymmetric 
contributions to both B(£i — > £jj) and Aa M (the simplified results in the mass-insertion 
approximations in Eqs. f|T5|) - f|T9|) have been shown only for the sake of clarity). The red 
areas are the regions where the 5-physics constraints are fulfilled. In our scenario the 
.B-physics constraints put a lower bound on M B and therefore, through the funnel-region 
relation, also on M 12 (see Figs. |3]and|I]). As a result, the allowed ranges for Aa /t and 
B[li — > Ejj) are correspondingly lowered. A complementary illustration of the interplay of 
B physics observables, dark- matter constraints, Aa M , and LFV rates -within our scenario- 
is shown in Figure ITlrl 

The normalization \5\ 2 L \ = 10 -4 used in Figures |6] and [7] corresponds to the central 
value in Eq. ffT4"j) for c u = 1 and M VR = 10 12 GeV. This normalization can be regarded 

9 For comparison, a detailed study of LFV transitions imposing dark-matter constraints -within the 
constrained MSSM with right-handed neutrinos- can be found in Ref. [45]. 
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Figure 7: Isolevel curves for B{fi — > ej) and £>(r — > /ry) assuming \5^ L \ = 10~ 4 and 
I^IlI = 10~ 2 i n ^e tan/3-M# plane. The green/red areas correspond to the allowed 
regions for the low-energy observables illustrated in Figure [3] for [pi, = [1.0,0.4] TeV 
(left plots), [/i,M z ] = [0.5,0.4] TeV (right plots). 
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Observable 


Exp. bound 


Bound on 
the eff coupl. 


Expected \5 LL \ in MLFV 
for M UR = 10 12 GeV 


B(n -> cy) 
B(t — > e 7) 


< 1.2 x l(T n 

< 1.1 x 1(T 7 

< 6.8 x 10" 8 


\S?L 

\8f L 


< 3 x 10~ 4 

< 8 x 1(T 2 

< 6 x 10~ 2 


(0.3-3) x 10" 4 
(0.3-3) x 10" 4 
0.8 x 10~ 3 



Table 1: Present experimental bounds on the radiative LFV decays of r and fi leptons 
[40] and corresponding bounds on the effective LFV couplings b l [ L . The bounds are 
obtained by means of Eq. (1191) setting Aa M = 20 x 10~ 10 . The expectations for the 8 l [ L 
reported in the last two columns correspond to MLFV ansatz in Eq. ()14p with c„ = 1 and 
M VR = 10 12 GeV. 

as a rather natural (or even pessimistic) choice. As can be seen from Figures [6] and [TJ 
for such natural choice of 5ll the fi — > ej branching ratio is in the 10~ 12 range, i.e. well 
within the reach of MEG [42] experiment. Note that this is a well-defined prediction of 
our scenario, where the connection between /1 — ► and Aa M allows us to substantially 
reduce the number of free parameters. In particular, the requirement of a supersymmetric 
contribution to Aa^ of O(10~ 9 ) forces a relatively light sparticle spectrum and moder- 
ate/large tan/5 values which both tend to enhance the LFV rates. This fact already allows 
to exclude values of 6}? L above 10~ 3 , for which £>(/i — > ej) would exceed the present exper- 
imental bound0 Within the MLFV hypothesis, this translates into a non-trivial upper 
bound on the right-handed neutrino mass: M UR < 10 13 GeV. 

On the other hand, the normalization \5]^ L \ = 10~ 2 adopted for the r — ► ^7 mode is 
more optimistic given the MLFV expectations in Table [H We have chosen this reference 
value because only for such large LFV entries the r — ► ^7 transition could be observed 
in the near future. From the comparison of Figure E] and Tabled] we deduce that, unless 
/i — > e7 is just below its present exclusion bound, an observation of r — > fij above 10~ 9 
would exclude the LFV pattern predicted by the MLFV hypothesis [37]. 



4 Conclusions 

Within the wide parameter space of the supersymmetric extensions of the SM, the regime 
of large tan j3 and heavy squarks represents an interesting corner. It is a region consistent 
with present data, where the (g — 2)^ anomaly and the upper bound on the Higgs boson 
mass could find a natural explanation. Moreover, this region could possibly be excluded 
or gain more credit with more precise data on a few 5-physics observables, such as 

10 For M VR 10 12 GeV other sources of LFV, such as the quark-induced terms in Grand Unified 
Theories cannot be neglected [41]. As a result, in many realistic scenarios it is not easy to suppress LFV 
entries in the slepton mass matrices below the 10 -4 level [38]. 

11 For a recent and detailed analysis on the bounds for LFV soft breaking term as functions of the 
relevant SUSY parameters (without assuming the present g — 2 anomaly as a hint of New Physics), see 
Ref. [46]. 
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BiB — > tv) and BiB — > £ + £~). In this paper we have analysed the correlations of the 
most interesting low-energy observables within this scenario, interpreting the (g — 2) M 
anomaly as the first hint of this scenario, and assuming that the relic density of a Bino- 
like LSP accommodates the observed dark matter distribution. In view of improved 
experimental searches of LFV decays, we have also analysed the expectations for the rare 
decays fi — > ej and r — > /i(e)7 in this framework. 

The main conclusions of our analysis can be summarised as follows: 

• Within this region it is quite natural to fulfill the dark-matter constraints thanks 
to the resonance enhancement of the X1X1 ~^ H,A — > // cross section (A-funnel 
region). As shown in Fig. [21 this mechanism is successful in a sufficiently wide area 
of the parameter space. 

• From the phenomenological point of view, the most significant impact of the dark- 
matter constraints is the non-trivial interplay between and the 5-physics ob- 
servables. A supersymmetric contribution to a M of O(10~ 9 ) is perfectly compatible 
with the present constraints from BiB — > X a j), especially for Ajj < 0. However, 
taking into account the correlation between neutralino and charged-Higgs masses 
occurring in the A-funnel region, this implies a sizable suppression of BiB — > tv) 
with respect to its SM prediction. As shown in Figure [5j the size of this suppression 
depends on the slepton mass, which in turn controll the size of the supersymmetric 
contribution to a^. In particular, we find that Aa^ > 2x 10 -9 implies a relative sup- 
pression of BiB — ► tv) larger than 10% A more precise determination of BiB — > tv) 
is therefore a key element to test this scenario. 

• A general feature of supersymmetric models is a strong correlation between Aa^ and 
the rate of the LFV transitions ii — > ij'j [34]. We have re-analysed this correlation 
in our framework, taking into account the updated constraints on Aa^ and B- 
physics observables, and employing the MLFV ansatz [37] to relate the flavour- 
violating entries in the slepton mass matrices to the observed neutrino mass matrix. 
According to the latter (pessimistic) hypothesis, we find that the /x — > branching 
ratio is likely to be within the reach of MEG [42] experiment, while LFV decays of 
the r leptons are unlikely to exceed the 10~ 9 level. 
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